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Within the framework of the QED3 theory of cuprate superconductivity it is argued that bulk 
antiferromagnetic (AF) order can coexist with d-wave superconductivity in the underdoped region, 
in agreement with recent experiments. The AF order arises from the phase fluctuating d-wave super- 
conductor via the mechanism of spontaneous chiral symmetry breaking, provided that fluctuations 
are sufficiently strong. The phase diagram for this coexistence is mapped out by means of analytical 
and numerical solutions of the underlying Dyson- Schwinger equation in the large N limit. 



Interplay between the AF order and superconductiv- 
ity remains one of the central themes in the physics of 
the high-T c cuprates 0, 0, 0, El- Recently, an appeal- 
ing new connection between AF and d-wave supercon- 
ducting (dSG) orders has been pioneered, based on the 
ideas originally articulated by Emery and Kivelson 
When long range dSC order is destroyed by thermal or 
quantum vortcx-antivortex fluctuations [S H H> 0> Ell j 
the resulting state can be either a symmetric algebraic 
Fermi liquid or, if the fluctuations are sufficie ntly 
strong, an incommensurate antiferromagnet [T2I IT^L ITlj 
with ordering vector Q illustrated in Fig. 2] In the lat- 
ter case AF order arises through an inherent dynamical 
instability of the underlying effective low energy theory 
of a phase fluctuating ci-wave superconductor, a (2+1) 
dimensional quantum electrodynamics, QED3 |lCt lll| . 
This instability is known as the spontaneous chiral sym- 
metry breaking (CSB) [lj, [lfl [lTj and is a well studied 
phenomenon in the particle physics. 

In this communication we investigate the possibility 
that bulk AF order can set in within the superconduct- 
ing phase, resulting in the region of coexistence of AF 
and dSC orders in the phase diagram shown in Fig. ^ 
within the QED3 framework. Other approaches to this 
problem have been considered previously; see e.g. Ref. 
|l8| . Recently, experiments have found tantalizing hints 
of such coexistence in zero app lied magnetic field in Y 
and La based cuprates 0, 120L l2ll I22I l23| . It has been 
shown previously [24| that within the QED3 framework 
such coexistence indeed can occur locally in the vicinity 
of fluctuating field-induced vortices, as found in neutron 
scattering ill HI, /zSR ^ and STM experiments. 

Within the QED3 theory [Tjj the dynamical agent re- 
sponsible for the emergence of AF order is a noncompact 
U(l) Berry gauge field which encodes the topological 
frustration encountered by the nodal fermions as they 
propagate on the background of the fluctuating vortex- 
antivortex plasma. In the non-superconducting phase 
Berry gauge field is massless 0, 0, 0|. Its quanta, 
"berryons" , have the same effect as photons in ordinary 
QED3: they mediate long range interactions between 
fermions and precipitate the chiral instability if the num- 
ber of fermion species N (pairs of Dirac nodes per unit 
cell in a dSC) is less than a critical value iV c ~ 3.2 [l5l[u^ |. 
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FIG. I: a) Schematic phase diagram of a cuprate super- 
conductor with AF/dSC coexistence. AFL denotes algebraic 
Fermi liquid, the symmetric phase of QED3. b) Cuprate 
Fermi surface with internodal wavevector Q. 



In the dSC phase, as vortices bind to finite pairs or loops, 
the Berry gauge field becomes massive [l(J, HH Q-3 • Nom- 
inaly, one would expect that coupling to such massive 
gauge field should become irrelevant for the low energy 
physics. Here we show that this is not necessarily so. 
Based on the approximate analytical and full numeri- 
cal solutions of the underlying large- N Dyson-Schwinger 
equation for the fermion self energy we find that solutions 
with broken chiral symmetry and finite fermion mass m d 
can be found even when the gauge field acquires small 
mass m a ■ This opens up a possibility for the coexistence 
of bulk AF and dSC orders in the QED3 theory of un- 
derdoped cuprates. 

As shown in Refs. jToL HH the low-energy dynamics 
of fermionic quasiparticles in a d-wave superconductor 
coupled to fluctuating vortices is governed by the QED 3 
Euclidean action S — J d 3 xCo with 

JV 

C D = *i{x)^{idu - O*i0*0 +C B [a{x)\. (1) 
1=1 

Here, ^i(x) is a four component Dirac spinor represent- 
ing the "topological" fermion excitations associated with 
a pair of antipodal nodes, x = (r, r) denotes the space- 
time coordinate, and 7^ (// = 0,1,2) are the gamma 
matrices satisfying {7^,7,,} = 2<5 M „. The number TV of 
fermion species is equal to 2nc u o, with nc u O denoting 
the number of CuO planes per unit cell ^lj . Lagrangian 
Lb encodes the dynamics of the gauge field a M and is 
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given by C B [a] = \U llv {q)a ll {q)a v (-q) with 



R»u(q) = (ml + a\q\ + q 2 ) 



(2) 



The mass and Maxwell terms in Eq. J3J) originate 
from the bare vortex action; m a — in the non- 
superconducting phase, reflecting the unbound nature 
of vortex-antivortex fluctuations, while m a > in the 
superconductor [ToL ITU H^ . The q- linear term, with 
a = Ne 2 /8, reflects the partial screening of the gauge 
field by the medium of topological fermions, expressed to 
one loop order. 

We analyze the CSB phenomenon in the Lagrangian 
(11121) by solving the Dyson-Schwinger (DS) equation for 
the fermion self energy E(p) to leading order in 1/N ex- 
pansion. This technique is known to be robust against 
the hig her order 1 /N corrections even for small N ~ 2 
0,123, and the results appear to be insensitive to vertex 
corrections |3fj|. 

With the above definitions we can write the DS equa- 
tion for the fermion self energy as 
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d 3 k 7m^(p ~ *)S(fc)7v 



(2tt) 3 



E 2 (fc) 



(3) 



where D fJlU (q) = IT" 1 (q) is the gauge field propagator 
in the Landau gauge. Following Pisarski [13, the sim- 
plest way to establish the spontaneous mass generation 
for fermions is to assume that, crudely, E(p) is constant 
over most of the integration range, E(p) = mp, and de- 
mand self consistency for the constant value ran'- 



mD= N 



a f d 3 k ~f ll D llv (k)mD')v 



(2n) 



(4) 



The above DS equation always has a trivial solution 
nip = 0, corresponding to the chirally symmetric state. 
A non-trivial, symmetry broken solution mo > satisfies 



thus suggests that CSB in QED3 could survive in the 
presence of small gauge field mass. Translated into the 
language of condensed matter physics Eq. © suggests 
that even bound vortex-antivortex fluctuations could give 
rise to AF instability within the state with the true dSC 
long range order. 

One can also perform the integral appearing in Eq. (J5J 
exactly and confirm that the solution ijfJJl indeed emerges 
in the limit m 2 <C amo (up to an unimportant prefac- 
tor). The approximation employed to obtain Eq. IjBJ is 
instructive in that it reveals the fundamental reason for 
the insensitivity of CSB to small gauge field mass: the 
largest contribution to the RHS of Eq. (JSJ comes from 
the region mo < k < a. Therefore, modifying the gauge 
boson propagator at momenta k <C mo by introducing 
the mass term has very little effect on the system. 

The above solution is highly suggestive but because 
of the crude approximations involved it is deficient in at 
least two ways. First, for m a = 0, it fails to reproduce the 
result of a more refined treatment of DS Eq. J3J ,3, 
that no CSB occurs for N > N c , with iV c = 32/tt 2 to 
leading order in 1/N ^(|. Second, Eq. iJBJ is only valid 
for m 2 <C amn- In what follows we present a more 
careful treatment of Eq. © which is free of the above 
deficiencies. We use these solutions to map out the phase 
diagram of CSB as a function of N and m a . 

To this end we must relax our assumption of constant 
E(p) in Eq. ((3J) and treat its full functional dependence 
on the three-momentum p. ft is still possible to perform 
the angular integral to obtain 
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E 2 (fc) 
k+p + 02 

\k-p\ + e 2 j\ 



where 61,2 — \{& ± \J oi 2 — 4m 2 ). In the limit a 2 ^> m 2 
we recover the equation studied in Ref. |16| . 



I = 



k 2 dk 
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tt 2 NJ (k 2 +m 2 D )(ml + ak + k' 2 )' 



(5) 



where we have taken the trace of both sides and per- 
formed the angular integrals. For small gauge field mass, 
m 2 <C arriD, we observe that the fermion mass effectively 
cuts off the integral in the infrared while a plays the role 
of the ultraviolet cutoff. We may thus approximate Eq. 
© by 



I = 
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k 2 dk 



■In 



n 2 N J mD k 2 (ak + mfj n 2 N \amD+m 2 
In this approximation the fermion mass is given by 



mo = ae 



*N/8 



m- 



(G) 



This is just the classic result of Pisarski [15| with a small 
correction due to the gauge field being massive. Eq. © 



ir 2 Np 



dk 



kT,(k) 
fc 2 + E 2 (fc) 



i(k,p), (8) 



but with modified lower bound on the integral. In the 
massless theory (m a = 0), Eq. © is solved by linearizing 
the integrand in E and seeking the solution in the form 
of a power law, E(p) ~ p b . Only solutions with complex 
b turn out to be physically admissible placing a condition 
on the number of fermion species N < N c = 32/tt 2 . The 
chiral mass is then given by [l6| 



m D = E(0) 



cae 



-2ir/yjN a /N-l 



(9) 



with c a numerical constant. 

Returning to the theory with massive gauge field, we 
observe that a simple power law ansatz no longer solves 
the linearized integral equation ©, essentially because 
of the appearance of the nonzero lower bound of the in- 
tegral. Consequently, it is a non-trivial issue to extend 
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FIG. 2: Numerical solution of Eq. (0 for dynamically gener- 
ated fermion mass mo as a function of gauge field mass m a 
for selected values of N < N c . Inset: curves for various N 
scaled as described in the text. Dashed line represents the 
solution to Eq. iftTty. 



such an analysis to our case of interest with m a > |3lj | . 
Fortunately, it is relatively straightforward to analyze the 
full nonlinear, angle integrated, DS Eq. (JJJ numerically. 
We start from a constant £(&) and numerically evaluate 
the RHS at n discrete points pi. We then iterate this 
procedure until the solution for S(p) no longer changes 
appreciably between iterations. In this we adopt an up- 
per cutoff, A, for the integral. We find that when A > a, 
our solutions no longer depend on A. 

Figures|21and|21summarize our findings. Fig.|21displays 
the dependence of chiral mass on the gauge held mass 
m a for selected values of N < N c . A notable feature here 
is the surprising resilience of the chiral mass: mjj persists 
until very large values of m a . A hint of why this may be 
so is offered by a simple power counting in Eq. JSJ. This 
reveals that the relevant scale actually is not m a itself but 
m^/a, with m a /a <C 1. The latter quantity also enters 
our estimate for iriu, Eq. © , which Fig. [21 confirms to 
be valid for small m a . 

When scaled by m* D and m*, defined as endpoints of 
the curves in Fig. El the points for different N fall on the 
same universal curve (see inset). One way to understand 
this universality is from Eq. iJ^J: its exact solution is given 
by the implicit equation 



\m D 



0, 



(10) 



with f(x,y) = x 2 hix + y 2 hiy + ^xy, independent of 
N. This solution is displayed as a dashed line in the 
inset of Fig. J2J). We find that the scaling holds very 
well up to N w 2.9; closer to AT° = N c (m a = 0) = 32/tt 2 
the agreement gets progressively worse. This comparison 
indicates that the approximation of constant self energy, 
leading to Eq. (|SJ), is very accurate for determining mo, 
as long as N is not too close to N®. 

Another way of viewing the above data is to consider 
N c as a function of m a . Fig. |3] shows this dependence 
as determined by our numerical solution of Eq. (JJJ. As 
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FIG. 3: Numerical solution of Eq. for N c as a function of 
gauge field mass m a . The dashed line represents a fit to the 
heuristic formula (II 12 P with C = 2.7. 



expected N c is a decreasing function of m a . A striking 
feature is that N c remains relatively large for m a up to 
significant fraction of a. The following simple argument 
leads to a highly accurate heuristic formula for N c (m a ) 
At low energies there are only two scales in the problem: 
mo and m a . Increasing m a leads to reduction of mo- 
lt is then natural to expect that vanishes when the 
two scales cross, i.e. m a P3 m* D . Here m* D is the fermion 
mass for given N but with massless gauge field, given 
by Eq. @. Given our previous insight that the relevant 
quantity is m 2 /a we take our criterion to read 



(m a /a) 2 = C(m* D /a), 



(11) 



with C is a numerical constant. Inserting m* D from Eq. 
(J5J) and solving for N we obtain 
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N c (m a ) = -j 



7r/ In 



Co, 



(12) 



where C = y/cC. As seen from the fit in Fig.[3]for C sa 2.7 
this leads to an excellent agreement with our numerical 
results. We note that the same expression, with C = 
e 2 S3 7.4, was obtained by Gusynin et al. |3l| from the 
linearized version of the DS equation with an infrared 
cutoff. 

Our results presented above firmly establish the ex- 
istence of nontrivial, chiral symmetry broken solution of 
the QED3 Dyson-Schwinger equation in the presence of a 
small gauge field mass. The question arises whether this 
represents a true ground state of the system, as is known 
to be the case for the massless-photon QED3 [lfj. We 
believe this remains true when m a > 0. Imagine turning 
on the gauge field mass from zero deep inside the symme- 
try broken phase, where there is large energy difference 
between the ground state and the "false vacuum" repre- 
sented by chirally symmetric phase. It is then clear that 
introduction of arbitrarily small gauge field mass cannot 
raise the energy of the ground state above that of the false 
vacuum: the chirally broken phase must remain a ground 
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state of the system for some range of values m a < m*. 
as indicated by our solutions © and Fig. [21 

Detailed nature of the criticality in the presence of 
gauge field mass appears to be a nontrivial issue. The 
transition in the ordinary massless-photon QED3 is 
thought to be of a special "conformal" variety, i.e. a con- 
tinuous transition of infinite order |33l l34j . This is due 
to the long range nature of interactions mediated by the 
massless gauge field. One would therefore expect that in 
the massive case the transition becomes a conventional 
second order 0, . 

Recently, powerful field theoretic duality arguments 
have been advanced [35j leading to a proposal for an up- 
per bound N® < |. This suggestion finds some support 
in numerical simulations of noncompact lattice QED3 
|36j | which found no decisive signal for chiral mass genera- 
tion at N = 2. If these results are correct this would sug- 
gest that the conventional state of the art analysis based 
on the Dyson-Schwinger equation © overestimates N® 
by more than a factor of 2. Our analysis, based on this 
same technique, would then also be quantitatively inac- 
curate. However, we expect its qualitative features to 
remain valid. 

If we stipulate that chiral symmetry breaking occurs 
in the massive-photon QED3, does this necessarily im- 
ply coexistence of the AF and dSG orders in cuprates? 
Answering this question involves additional subtlety that 
is related to the detailed nature of the criticality at the 
transition to the dSC state. According to the standard 
phenomenology 0, 0, 0, 0] , approaching the transi- 
tion from the above (i.e. from the symmetric pseudo- 



gap phase) berryon is massless but the charge e tends to 
zero, transforming eventually the Maxwell term into the 
mass term in the dSG phase. Since according to Eq. © 
mr> is proportional to e 2 (through a) this would suggest 
that chiral mass should vanish at the transition. How- 
ever, there is presumably higher order (~ q A ) term in 
the bare berryon action whose prefactor also diverges at 
the transition giving rise to a Maxwell term in the dSC 
phase. The chiral mass would then be proportional to 
this new "charge" in the dSC phase. Our calculations 
above strongly suggest that such a chirally broken phase 
can be a globally stable ground state of a system with 
massive berryons. Again, while this work might have 
resolved the nature of the phases on both sides of the 
transition, the detailed understanding of the criticality 
itself appears to be a complex problem awaiting future 
resolution. 

QED3 theory as formulated in Refs. 0, 0] therefore 
appears to support a locus of AF/dSC coexistence in the 
low doping region of the phase diagram Fig. ^ In this re- 
gion the system remains superconducting while fcrmionic 
excitations become fully gapped. This small gap should 
be observable in thermodynamic and transport measure- 
ments. In particular the superfluid density should be- 
come exponentially activated in very clean samples at 
low temperatures. 

The authors are indebted to I. F. Herbut, D. E. Sheehy 
and Z. Tesanovic for many penetrating discussions. This 
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